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On the Solutions of Certain Types of Linear Differential 
Equations with Periodic Coefficients. 

By F. R Moulton* and W. D. MacMillan. 



§ 1. Introduction. 

Most of the methods which are employed for finding the solutions of differ- 
ential equations were devised in order to solve the practical problems which 
arise in celestial mechanics. It is sufficient to mention in this connection the 
expansion of the solutions as power series in the independent variable, as power 
series in parameters, the method of the variation of parameters in the general 
non-linear case, and the method of successive approximations. The first two 
and the last were used formally by the founders of the analytic theory of the 
motions of the planets — Clairaut, d'Alembert, and Euler — and the third was 
given its general formulation by Lagrange at the end of the eighteenth century, 
and its widest application by Delaunay in his Lunar Theory, in the middle of 
the nineteenth century. All of these processes were extensively employed in 
celestial mechanics for obtaining practical results, without any inquiry being 
made regarding the circumstances and realm of their validity. Indeed, it was 
as late as 1842 that Cauchyf began laying the foundations of the modern theories 
of differential equations. 

The recent contributions to the theory of differential equations, which have 
been stimulated by problems in celestial mechanics, have come chiefly from the 
hands of Hill and Poincare. In 1877 Hill J published privately at Cambridge, 
Mass., his famous investigation of the motion of the lunar perigee. In this 
memoir he treated with rare skill the differential equation 

^ + 6(t) w = o } 

* Research Associate of the Carnegie Institution of Washington. 

tCauchy's Collected Works, 3d series, Vol. VII. 

t Also reprinted in Acta Mathematica, Vol. VIII (1886), pp. 1-36; Hill's Collected Works, Vol. I, pp. 343-270. 
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where 6 is a simply periodic function of t. About the same time Hermite* dis- 
covered the form of the solution of Lamp's equation, which has a doubly periodic 
coefficient. Starting from Hermite's results Picard f showed that in general a 
fundamental set of solutions of a linear differential equation of the n-th order 
having doubly periodic coefficients of the first kind can be expressed in terms of 
doubly periodic functions of the second kind. In 1883 FloquetJ published a 
complete discussion of the character of the solutions of a homogeneous linear 
differential equation of the n-th order having simply periodic coefficients. In 
this memoir Floquet gave not only the form of the solution in general, but he 
considered in detail the forms of the solutions when the fundamental equation 
has multiple roots. The forms of the solutions being thus known, the efforts of 
later writers have been directed toward the discovery of practical means for 
their actual construction. Among those who have discussed the problem of 
finding the solutions of Hill's equation we may mention Lindemann §, Lind- 
stedt||, Bruns^J, Callandreau ** , Stieltjesff, and HarzerJJ. 

In all of these investigations a large amount of attention has been devoted 
to finding the roots of the fundamental equation, or equivalent transcendentals. 
Hill determined them from an infinite determinant which he first introduced 
into analysis in this connection ; Lindstedt found them from an infinite con- 
tinued fraction §§. In all cases these transcendentals were computed first, and 
then the solutions were found later. It should be noted also that the processes 
are valid only under certain special conditions which, fortunately, are satisfied 
in the case of Hill's equation. The problem is treated in a much more general 
way in Poincare's Les Methodes Nouvelles de la Mecanique CUeste, Vol. I, Chapter 
II, Sec. 29, and Chapter IV. 

In dynamical problems involving accelerations simultaneous differential 
equations of the second order naturally arise, but it is easy to reduce them to 
twice the number of simultaneous equations of the first order. Now n simul- 

* Comptes Bendus, 1817 et seq. 

f Oomtes Bendus, 1879-1880 ; Journal fur Mathematik, Vol. XC (1881). 

J Annates de VEcole Normale SupSrieure, 1883-1884. 

§ Maihematisehe Annalen, Vol. XXII (1883), p. 117. 

|| Astronomische Nachrichten, No. 2503 (1883), and MSmoires de VAcadimie de St. PStersbourg, Vol. XXI, No. 4. 

If Astronomische Nachrichten, Nos. 2533 and 2553 (1883). 

**Ibid., No. 2547(1883). 

t+ Ibid., Nos. 2601 and 2609 (1884). 

tt Ibid., Nos. 2850 and 2851 (1888). 

§§ See Tisserand's Mecanique Celeste, Vol. Ill, Chapter I. 
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taneous differential equations of the first order include one differential equation 
of the n-th order, for the latter can always be reduced to the former. 
Hence we shall treat here as the general case and the one most simply 
connecting with dynamical problems a set of simultaneous linear differential 
equations of the first order having simply periodic coefficients. We shall find 
the character of the solutions of the differential equations without further 
restrictions by a very direct process. Then, simple and convenient methods 
are given for constructing the solutions in all cases in which the coefficients of 
the differential equations are expansible as power series in a parameter p, and 
the terms not depending upon (i (at least in a large part of the discussion) are 
constants with respect to the independent variable. The linear differential 
equations having periodic coefficients which arise in celestial mechanics,* of 
which Hill's equation is a simple example, belong to this class. 

§ 2. The Fundamental Equation. 
We shall consider the differential equations 

n 

x[ = 2 dtj(t)xj, i—\, , n, (1 ) 

where x[ is the derivative of x t with respect to the independent variable t, and 
where the 0^ are uniform analytic functions of t and are periodic with the 
period 2n. Let 

*ii = $a(t), , x in = <M0> * = 1, ,n, 

be a fundamental set of solutions of (l), where x v = ^) v (t), i= 1, . . . . ,n, is the 
j-ih solution. The determinant of the fundamental set, 

A = |<M, 

is found, by taking the derivative and reducing by means of (1), to satisfy the 

relation „ 

A' = A2 0„(«); 

whence f 

. 2 en at 
A = A>e °"' (2) 

Hence A can become zero or infinite only at a singularity of some 6 u (t). 

*Some of the methods exhibited here were devised in connection with problems raised in the theory of 
periodic orbits. They were first applied to Hill's equation by Moulton in a paper whose abstract is in Bull, 
of the Am. Math. Soc, Vol. XIII (1906-7), p. 71, and were later extended by our joint investigations to the 
general case. 

fDarboux, Comptes Rendus, Vol, XC (1880), p. 536, 

9 
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We now start from a particular set of solutions satisfying the initial con- 
ditions 

* M (o) = i, <Mo) = o, if y 4= *• ( 3 ) 

It is clear that a set of n solutions satisfying these initial conditions can be con- 
structed, and since for them A = 1 is distinct from zero, they form a, fundamental 
set of solutions. 

Let us make the transformation 

x t = e a %, (4) 

a being an undetermined constant. Then equations (1) become 

n 

yi + ay i = 2 Bijyj. (5) 

Any solution of (5) can be written in the form 

n 

y t = e-* 2 A^it), i = 1, . . . . , n, (6) 

where the A$ are suitably chosen constants. 

We now inquire whether it is possible to determine a and the Aj so that 
the y i} as defined by (6), shall be periodic with the period 2n. From the form 
of (5) it is clear that sufficient conditions for the periodicity of the y t with the 
period 2n are 

y,(270-2fc(0) = O, i=l,....,n. (7) 

Imposing these conditions on (6), we get 

2 4[fc/2») - e^/O)] = 0, • = 1, . . . . , n. (8) 

In order that these equations may have a solution other than J. x = .... = A n =0, 
the determinant of the coefficients of the Aj must equal zero. Making use of (3), 
representing 4> M (27t) simply by <p ij} and letting e 2a,r = s, the determinant is 

#11 — s > #12 , , $m 

$21 ) $22 s ) • • • • , <p2n 



D = 



> > • • • •> • • • • 

tyn\ > tynZ > • • ■ * j <Pnn s 



0. (9) 



This is the fundamental equation associated with the period 2n and is of degree n 
in s. It admits neither s = nor s = oo as a root, since the absolute term is the 
value of the determinant of the fundamental set of solutions at the regular point 
t = 27i, and the coefficient of s n is ( — 1)". 
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§ 3. Form of the Solutions. 

Suppose the roots of (9) are s 1 , . . . . , s n and that they are all distinct. For 
each of them there is therefore at least one first minor of D which is distinct from 
zero, and hence for each of them the ratios of the Aj can be determined from (8) 
so that the y i shall be periodic. In this way n solutions are obtained which can 
be shown to constitute a fundamental set. 

Suppose now that s 2 = s t and that all the other roots of (9) are distinct. 
There are two cases to be considered, according as all of the first minors of D 
vanish or do not vanish for s = s x . Suppose all the first minors of D vanish 
for 8=$!; since by hypothesis s = Sj is a double root and not a triple root, 
there is at least one second minor of D which does not vanish for s = s t . 
Therefore two of the A$ can be taken arbitrarily and the remaining n — 2 can 
be expressed in terms of them so that the y % shall be periodic. We thus obtain 
two distinct solutions of the form 

x a = e a *y a and x i2 = e«*y i% , i — 1, . . . . , n, (10) 

where the y a and y i2 are periodic. 

Suppose s = s 1 is a double root of (9) and that not all its first minors vanish 
for 8 = 8!. Then there is but one solution of the form 

x a = e^y a . 

The y a are expressible linearly in terms of the $ i3 by (6). Let the notation be 
chosen so that the minor which is not zero is formed from the elements of the 
last n — 1 columns. Then A x must be distinct from zero in order to avoid the 
trivial case in which all the A 5 are zero. 
Then we take as a new set of solutions 

<% = e a *y iU x ij = ^) ij (t), i=l, , n;j=2, , n. (11) 

These solutions constitute a fundamental set, for their determinant is 



e ai t 



I/u, <?>i2 > > $m 

V%\, Qnz, , <frs» 

. . . , . . ., . . . ., ... 

Vn\, $n2, , Qnn 



which becomes, by means of (6), .AJ^yl, which is distinct from zero. Starting 
with this fundamental set, the fundamental equation is 
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D = (s-s 1 )D 1 = (s-s 1 ) 



2/ii(0) , $,„ , 



<?>ln 
<?>2ra 



= 0. 



(12) 



2/»l(°)> #»2 , , $nn — S 

Since s= St is a double root of D, the determinant Z^ has a single factor s — s v 
Since (11) constitute a fundamental set, any solution can be expressed in the form 



xt = Byfr*y a + X Bjfyj, 



i=l. 



j=3 



• 1 



n. 



(13) 



Now we make the transformation, corresponding to (4), to get a second 
solution associated with a lf 

a* = «**(y tt + ty fl ). (14) 

Imposing the condition that the x i% shall satisfy (l) we find, since e ait y a is a 
solution, 

n 

2/4 + a i tfa = 2 0« y i% — y iU i = 1, . 



j=i 



. , n. 



Therefore sufficient conditions that the y® shall be periodic with the period In 
are 

n 

y i2 (2n) — y iS (0) = = — 2*y fl (0) + 2 B } \f+*frtfn) - ^0)]. 

.7=2 

Substituting Sj for e 2 " 1 '', we have 

- 2*9,3^(0) + 2 4[4>«(2t) — «i^(0)] =0, » = 1, . 

i=2 

The condition that these equations shall be consistent is 



, n. 



(15) 



A = 



2/11(0) , 

2/21(0) , 



#12 

#22 " 



'1) 



2/nl(0), fc 



'm2 



4»m 

<?>2n 



= 0. 



It was shown in (12) that this determinant has a single root s = s x . Hence not 
all its first minors are zero. By hypothesis not all minors of the first order 
formed from the last n — 1 columns of D, which are the same as the last n — 1 
columns of D lf are zero. Therefore we can solve equations (15) for B % , . . . . , B n 
in terms of y a (0) which carry an arbitrary constant as a factor. Consequently 
in this case a second solution depending upon a x exists and is expressible in the 
form (14), where y a and y i% are both periodic. 
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When s = Sj is a triple root of D = 0, an analogous discussion shows that 
the three associated solutions are of the form 



or 



or 



*n = ^%i, *® = ^*y«, «« = *-*[>« + %*i + y«)]> t = 1, . . . . , n, 
«a = « a, Vii> 35« = ^Oe + fy*J, «« = e a ' e [y i3 + ty® + l t*y a ~\, i = 1, , n, 



according as all the minors of D of the first and second orders vanish for s = s 1} 
or all those of the first order but not all those of the second order vanish, or not 
all the minors of the first order vanish. 

In case s = s t is a root of multiplicity v, the associated group of solutions 
have the form 

x il = e^ t y il , 

a« = **[&« + ty* + l <ty«] , (1 6) 



x. 



= <**[[>„ + ty^ + ....+ ( V _L 1){ f-'yn]. 



where the y i} are periodic with the period In, provided not all the minors of D 
of the first order vanish for s = s 1 . In case all the minors of D of order 1c — 1, 
but not all of order h, vanish for a = 8 lt there are k solutions of the type of 
the first of (16), the others involving products of t to powers not exceeding 

v — Tc and the y tj , *=1, , n, j = 1, . ; . ., v. The details for a single 

differential equation of order n were given by Floquet (foe. cit.), and the results 
are similar here. 

§ 4. The Characteristic Equation when the Coefficients, i;j , are Power Series 

in a Parameter (i. 

We shall now assume that the 6 %i are expansible as power series in a 
parameter (i whose coefficients are separately periodic with the period 2n, and 
that the series converge for all real finite values of t if | (i | < p. We shall assume 
further that for (i = 0, W = a i} , where the a tj are constants. Under these con- 
ditions, which are often realized in practice, particularly in celestial mechanics,* 

* Hill's differential equation, considered first in his paper on the motion of the lunar perigee, is an 
equation of this type. 
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the discussion of the character of the solutions can be made so as to lead to a 
convenient method for their explicit construction. 
Consider the equations 



n n oo 

as,' = 2 0« xj = 2 [« y + 2 flgyfo , 

^=1 5=1 A=l 



• = 1, . ..., n, 



(17) 



where the a i7 - are constants. When p = they have the solution 






where a (0) is any one of the roots of the characteristic equation 



— a<°> 



a. 



31 



a 



n\ 



«12 



'■13,' 



a 



(0) 



a 



n% 



• * J 

• ' i 



*"Zn 



._«.«> 



= 0. 



(18) 



If the rc roots of this equation are distinct, a fundamental set of solutions of 

(17) for ^ = is 

3.(0) — r ffl.m j ,• — 1 „ 

wnere the c i} can be taken so that \c i} \ = 1. 

By Poincare's extension of Cauchy's theorem,* equations (17) can be 
integrated as power series in ft which will converge for < t < T, T any arbi- 
trary finite value, provided \(i\<^E, E depending upon T. Hence we can write 
the general solution of (17) in the form 



x t = 2 ■Afci/rf* + 2 ztftt*], 



(19) 



where the A s are the constants of integration, and the x'ff are functions of t 
depending on the Q i} . We can take the initial conditions so that 



Therefore 



a5y(0) = 2a$>(0)p* = fly. 

k=0 



°®(0) = o ti , xfl(0) = 0, 
Now we make the transformation 



h— 1, . 



»« = e at ^, 



after which the differential equations and their solutions become 



*Les MHhotXes Nouvelles de la M&canique Oileste, Vol. I, Chapter II. 
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y[ + a y t = 2 [«y + 2 0gy]2/«, » = 1, , n ; 

i=i &=i 

n oo 



(20) 



Imposing the conditions that the y t shall be periodic with the period 2n, viz., 
y i (2n)—y i (0) = Q, we get 

= 2 A£c t £M n -*'— 1) + e-*" 2 xf^n)^, i = 1, . . . . , n. (21) 

In order not to have the trivial case where all the Aj are zero, we must set the 
determinant 

D = | [e w (eW- ) " — 1) + <r»" 2 asg>(2«y] | = 0, (22) 

k=l 

which is a condition upon the undetermined constant a. This equation has an 
infinite number of solutions ; for if it is satisfied by a = a , it is also satisfied by 
a = a + vnV — 1, where v is any integer. All distinct solutions can be ob- 
tained with v equal to zero ; the others amount simply to taking periodic factors 
from the y i . The fundamental equation corresponding to (22) is obtained by 
putting e 2 "* = s. If the n values of s satisfying the fundamental equation are 
distinct, the corresponding values of a are distinct, but not necessarily the con- 
verse. "We shall use only those n values of a which reduce to the aj 0) for fi = 0, 
the n aj 0) being uniquely determined by (18). 

In the event that two of the roots of the characteristic equation, say af } 
and af\ are equal, the solutions are either of the form (19), where c4 0) = c4 0) , or 

Xi = A x [c a e^ + 2 x at f] + A[(c i2 + tc a )e^ + 2 «4V] 

fc=l k=l 

+ 2 Aifa/rf* + 2 a$y]. (23) 

3=3 *=1 

After making the transformation x i = e^y^ writing out the solutions for y t from 
(23), and imposing the periodicity conditions, y^n) — 2^(0) = 0, on the y i , we 
find that if not all the A } are zero either (22) is satisfied or 

D = [ [c it (e 3 H >-^_ l) + *+» | a@»(2ji)^] > 

[o e (^ a f- a) " — 1) + 27tc il e^i 0, -» ) " + e- 2a " 2 x%\2n) i i k ~\, . . . . | = 0, (24) 

k=\ 
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according as the solutions for ^ = are of the form (19) or (23), where the 
terms of the determinant not written are of the same form as those in (22). 
When for fi = the characteristic equation has a root of higher order of multi- 
plicity, a corresponding discussion must be made. If for ^ = the characteristic 
equation has several multiple roots, a corresponding discussion must be made 
for each associated group of solutions. 

§5. Solutions when af* Are Distinct and af } — aj 0) g£0 mod*/ — 1. 
The part of (22) independent of ^ is 

n 

D =\c i j(e*«P-^-l)\ = \c ij \n{e*«r-*»-l), 

j=i 

and the initial conditions have been taken so that |c w |=l. If (22) were an 
identity in (i, its n solutions would be the n solutions of D = 0, viz., a = aj 0) . 
In the general case in which it is not an identity in [t let 

a = 4 0) + /?'*• (25) 

Then we get 

n 

d = d, + pF^, ii) = ( e -^- -i)n'( e ^i 0) -r-^> _ i) 

+ [iF k (P k , fi) = 0, j$l; (26) 

where F k (/3 k , (i) is a power series in fj. and (3 k , converging for 

|&l<«, |pl<P><>- 
Since by hypothesis no aj 0) — a ( k 0) equals an imaginary integer, the expansion of 
(26) as a power series in @ k and (i contains a term in ft k of the first degree and 
no term independent of both fi k and (i. Therefore by the theory of implicit 
functions (26) can be solved uniquely for (3 k as a power series in ^ of the form 

8 k = (iP k (ix), (27) 

which converges for |^|>0 but sufficiently small. When we substitute this 
value of a = o4 0) + @k back in (21), we have n homogeneous linear equations 
among the A } whose determinant is zero, but whose first minors, for fi = 0, 
are not all zero ; since, by hypothesis, for p. = the roots of the determinant 
set equal to zero are all distinct and no two differ by an imaginary integer. 
Therefore the ratios of the Ay are uniquely determined by these equations as 
power series in (i, converging for | fi | sufficiently small. Substituting the ratios 
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of the Aj in (20), we have the particular solution y ik , » = 1, . . . ., n, carrying one 
arbitrary constant and expanded as a power series in p. Hence we may write it 

y« = «»2^. (28) 

Since the periodicity conditions have been satisfied, 

yjt + tot) - y ik {t) = 2 \tfi{t + tot) - y<i{tW = 

.7=0 

for every p whose modulus is sufficiently small and for all t. Therefore 

yg\t + tot) - yO](t) = 0, j = 0, ....,«, , 

from which it follows that each y\§, j = 0, . . . . , oo , is separately periodic. A 
solution is found similarly for each aj 0) . 

§6. Solutions when no Two af> Are Equal but when af — a^ 0) =0 mod V — 1 

Suppose two roots of the characteristic equation (18) for p = 0, say af and 
af\ differ by an imaginary integer and that there is no other such congruence 
among the aj 0) . Then the equation corresponding to (26) becomes 

n 

D = (e-^ - If II (eW-F-A). _ 1) + fr^(ft, p) + p 2 F^ u p) = 0. (29) 

J'=3 

The term of lowest degree in /3 a alone is +47t 2 /3f. The terms independent of (3 
carry p z as a factor, for every element of the first two columns of D, equation 
(24), in this case carries either /3 X or p as a factor. In order to get the terms in 
p alone we suppress those involving @ lf after which we get a factor p from each 
of the first two columns. In general the term of lowest degree in p alone will 
be in this case of the second degree. 

In a similar manner if p of the af are congruent to zero mod V — 1, then 
the term of lowest degree in (3 t alone is of degree^, and in p alone it is at least 
of degree p. 

The problem of the form of the solution of (29) is one of implicit functions. 
Writing the first terms explicitly, we have 

@1 + Yn@\P + 7oa^ 2 + terms of higher degree = 0, 

where y n , y 02 , .... are constants independent of /3j and p. The quadratic 
terms may be factored and we get 

(A — bi(*)(Pi — hi 1 ) + terms of higher degree = 0. 
10 
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If b x and b 2 are distinct, the two solutions have the form 

Ai = he + PPM, 1 
p a = hii + i*PM, J 

where Pj and P 2 are power series in (i which converge if | (i | is sufficiently small. 
In this case the solutions are found precisely as in § 5. 

If J I = 6 3 , the character of the solution depends upon terms of higher 
degree. It will proceed in powers of ±. V/x. unless the Q<$ satisfy special con- 
ditions. But for special values of the 0$ it may proceed according to integral 
powers of p. We shall consider in detail the case where the series is in powers 
of ±V(i. 

We see that the expansion of a x as a power series in Vfi will contain no 
term in V/x to the first power, but will have the form 

aj = af + Op* + afty + aJV + . . . . 

Suppose this series has been obtained from equation (29). Then since not all of 
the first minors of D are zero the ratios of the Aj will be determined from (21). 
Suppose [i8 is a non-vanishing first minor of D formed from the elements of its 
last n — 1 columns. Then it follows from the form of (22) that, solving (21), 
we get 

. __ fiS 2 A . (i 2 8j . . 



where 



^•r + sjv + ajV + ...-, y=2, 



Substituting these series for the A s in (20), we find that the y tl are developable as 
series of the form 

ya=Va + ygy + ySV + — , * = 1, — , «• 

Therefore the y n carry terms in ^ 4 , although the term in V/x is absent in the 
expression for o^ . 

If all the first minors obtained from D when the first column is suppressed 
are zero, and if there is a first minor which is distinct from zero when the second 
column is suppressed, the results are precisely the same. But suppose all the 
first minors obtained by omitting in turn the first and second columns are zero. 
Since the determinant has simple roots, there is at least one minor of the first 
order which is distinct from zero. Suppose it is obtained when the &th column 
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is suppressed. It follows from (22) that when a =-a x it will carry the factor (i z ; 
let it be [/8. Then solving (21), we gel 

a—HKa a — i*k a a— A a »■-* „ 

where 8 lf $ 2 , 8j do not in general vanish for ^ = 0. It follows from the first 
two equations that A k must carry ft as a factor. Hence in this case the y iX have 
the same form as before. The y iz have the same properties. These properties 
of the y a and the y i2 are necessary for the construction of the solution. 

We now return to the consideration of (29). If the discriminant of the 
quadratic terms of (29), 

D = r!i — 4^, 

is distinct from zero, the solutions are in integral powers of fz, and at least one 
of them starts with a term of the first degree in fi. But if the discriminant is 
zero, then the character of the solutions depends upon the coefficients of terms 
of higher degree. They may be either in integral powers of n or in powers 
of \//i . If the solutions are in V(i, they are real when (i has one sign and 
complex when it has the other. But if the solutions are in integral powers of n, 
they are either real or complex for both positive and negative values of (t. 

In all cases we get two solutions associated with the root a{ 0) . We should 
also get two solutions if we started from the root af\ but it follows from the 
form of (22) that they would not differ from those obtained by starting from af 1 . 
For each aj 0) + & the ratios of the A, are determined from (21), and the results 
substituted in (20) give the y^ The solutions associated with c4 0) , . . . ., a^' are 
found as in the preceding case. If there are several groups of a (0) in which these 
congruences exist, the discussion is made similarly for each one of them. 

§ 7. Solutions when <4 0) is a Multiple Root, 

Suppose that two and only two of the af\ viz. af and o4 0) , are equal, and 
that there are none of the congruences treated in §6. Then for ft = we get 
from (22) and (24) either 

A = ku(e K ^ _ " ) '-l), ObO^-^-I), ...., ^(e^f-^-l), ....|, 
or 

D = | c a {<?^-°* - 1), c i2 ( e s H°>-«>* _ i) + 27tc a e 2 <«f ->*, . . . . , 
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both of which, by the theory of determinants, reduce to 



n 



D = | c £ , | (e««f-)- - l) s n (<W->* - 1), ( | cy | = 1). (31) 

1=3 

If we let a = a£ 0) + A, & s before, and expand as a power series in (3 lt we 
find that the term of lowest degree in fa alone is 4n z fa. "When the determinant 
D is of the form (22) with o4 0) = a| 0) , the term of lowest degree in (i alone is at 
least of the second degree ; but when D is of the form (24), which is the general 
case, the term of lowest degree in ^ alone is in general of the first degree. 
Except in the special cases, the solutions for fa are therefore of the form 

A.= W), 1 

where P is a power series in p* containing a term independent of (i. The treat- 
ment of the special cases proceeds as in § 6. Substituting these expansions for 
a t = aj 0) + /?! in (21) the ratios of the A } are determined as power series in V/w, 
and these results substituted in (21) give the y a and y i2i as power series in V/t*. 

If, for (i = 0, p roots, ai 0) , , a£", are equal, then for these roots the 

expansion of D starts with fa as the term of lowest degree in fa alone, and, 
except in special cases corresponding to those above mentioned when two roots 
are equal, the term of lowest degree in (i alone is of the first degree. Conse- 
quently in general for af = . . . . = a ( ® we have 

faj = *V P ( e V), J = 0, • • • • , P — 1, 

where g is a ^>th root of unity. 

§8. Solutions when there Are Equalities and Congruences among the Boots of 

the Characteristic Equation. 

Suppose, for example, that a£ 0) = <4 0) and af differs from af> and of by an 
imaginary integer, and that there are no other equalities or congruences among 
the af\ There are two cases, (a) where the solutions are of the form (19) with 
af } = a4 0) , and (b) where the solutions are of the form (23). 

Case (a). In this case Z> '= becomes 

D = | c fl (^«P-»" - 1), % ( e *«?>-«>»_ 1), 0o (^-f>->- - 1), .... , 

e„(eW->- — 1), ....1 = 0. 
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Putting a = o4 0) + fa , we find that the term of lowest degree in fa is — 8n s fa. 
In order to get the term of Z> of lowest degree in (i we put /?j = 0. Then the 
first three columns of D are divisible by p, while the others do not contain ft as 
factor. Consequently the term of lowest degree in [i is of the third degree at 
least. Moreover, since the first three columns of D vanish with fa = (i = 0, 
there are no terms lower than the third degree in fa and /«. Hence in general 
D, in case (a), is of the form 

d = fa + y a /Sfc + n 2 /V + w 8 + •••• = o. (33) 

Since the coefficients of this equation are real, it always has at least one real 
solution for fa, vanishing with /i. The details of the various special cases are 
simply those of implicit functions. In general the three values of fa are ex- 
pansible in integral powers of p. 

Case (b). In this case D = becomes 

Z> = | Cfl(6*"F->"— 1), e («*P->*— 1) + 2nc a eW»-*' } c^eW-*- 1), . . . . | = o. 

By the theory of determinants this equation reduces to 

Z> = (e2(«r-«>* _ i)»( (S ««f»-")- — 1) n (e**P- a) — 1) = 0. 

j=i 

Introducing fa as before, we find that the term of lowest degree in fa alone is 
— 8it?fa. But when the terms involving (i are retained in D, the terms 
27tc a e 2(< f )-a) ' r can not be eliminated from the second column. Hence only the 
first two columns vanish for fa = ft = 0, and therefore in general in this case 
the expansion of D will contain a term in ft z alone. Since the first two columns 
vanish for fa = (i = 0, there will be no terms of degree lower than the second 
in fa and {i. Hence in general D, in case (b), has the form 

D = fa + ynfafi + yutf + .... = 0. (34) 

In the general case in which y n :£ and y 02 :£ there is one real solution in 
integral powers of ^ and two solutions in Vfi. The two latter are real or 
imaginary for p > according as y n is negative or positive. In all cases there 
is at least one real solution. 

The treatment of cases where there is a higher order of multiplicity of the 
roots <xj 0) and more numerous congruences is similar. If the total number of 
roots equal to aj[ 0) is v x , and of those congruent to a£ 0) mod s/ — 1 is v % , then in 
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the expansion of D the term of lowest degree in /3 a alone is of degree v x + v z , 
and the term of lowest degree in p alone is in general of degree v x . There are 
no terms of degree lower than v x in /?j and (i. 

§ 9. Solutions when D — Eas Two Boots Equal Identically in (i. 

The conditions necessary that D = shall have two or more roots identical 
in {i are that 

D(a )( i) = 0, JJL D (a,ii)=0 

for all \(i\ sufficiently small. Let us suppose, for simplicity, that a x and a 2 are 
identically equal in fi and that all the other a } are distinct. In this case D has 
the form (31), but there are no terms in D of the first degree in n, for then (3 n 
and p^ could not be equal. The value of @ n = /3 12 is found as in § 7, and the 
corresponding solution is obtained by solving (21) for the ratios of the Aj and 
substituting the results in (20). If all the first minors of D vanish for a = a!, 
then two of the Aj remain arbitrary and we obtain in this way the two solutions 
belonging to a x . 

Suppose not all the first minors of D are zero for a = a 1 . Then there 
remains but one arbitrary in the solution of (21) for the ratios of the Aj. Sub- 
stituting this result in (20) we get the solution y a ([i, t), from which we have by 
the general transformation 

^n = e ait y a , (35) 

where the y a are power series in fx. whose coefficients are periodic in t with the 
period 2n. 

By the general theory of § 2 we know that the other solution depending 
upon a x is 

*i3 = ^{y w + ty a ), i — 1, . . . • , n, (36) 

where the y i2 are periodic in t with the period 2n. Substituting these expres- 
sions in the differential equations (17) and making use of the fact that the e ttl *y a 
are a solution, we get 

n oo 

Va + aiVa ~ 2 [a tj + X CV*].%a = — Va, » = 1, , «■ (37) 

,?•=! k=l 

If the right members of these equations are put equal to zero, they become 
precisely of the form of the equations satisfied by y a . Consequently the only 
solution of these equations which is periodic with the period 2n is y i% = y a plus 
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such particular integrals that (37) shall be satisfied when the right members are 
retained. The part y a is useless since it belongs to the first solution, which con- 
tains an arbitrary factor. The method of finding the particular integrals will be 
taken up in § 16. It will there be shown that the y u are also power series in p. 
When D = has a multiple root of higher order of multiplicity for all \[i\ 
sufficiently small, the y Ui y i2 , . . . . , i = 1, . . . ., n, are found in succession, the 
problem for y i2 , y & , .... being that of linear equations with right members 
analogous to (37). 

§ 10. Direct Construction of the Solutions when the aS 0) Are Distinct and 
af — a| 0) g£ mod V— 1, » = 1, , n, j = l, , n. 

The methods given above lead to the constructions of the solutions, but the 
work is very laborious. However, the knowledge of their properties which we 
have obtained and their expansibility as power series in y. lead us to convenient 
methods for their construction. 

Under the hypotheses of this article it has been shown that there are n dis- 
tinct values of a, viz. a lt . . . ., a n , expansible as power series in (i such that 

Xik = ^ kt yik, i=l, , n, 

constitute a fundamental set of solutions, where the y ik are purely periodic and 
expansible as power series in p. Since the y ik are expansible as converging 
power series in /z and are periodic with the period 271, each y$ is separately 
periodic. Therefore the conditions on a = aj 0) + affy + • • • • and the y$ are 
that they shall satisfy the differential equations identically in p and 

S$(« + 2*)-j$(0 = 0. (38) 

The differential equations for the y t are 

n a, 

Vi + ay i = Z \a i} + 2 flJV] Vs* t = 1, . . . . , n. (39) 

For (i = the roots of the characteristic equation are a| 0) , . . . . , a<f . Consider 
any one of them, as a4 0) . Then a and the y t are expressible by converging power 
series of the form 



a = af + g4> + .... = X 4V » 



k=0 



Viu = tffi + V> + ••■■ = 2 2/&y, i = 1, 



., n. 

*=o 



(40) 
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Substituting (40) in (39) and equating coefficients of corresponding powers of [i, 
we have a series of sets of differential equations from which a k and the y ik can 
be determined so that the y ik shall be periodic with the period 2n. The deter- 
mination is unique except for the arbitrary constant factor of the solution. For 
simplicity this arbitrary will be determined so that yyj(0) = <hk' If c u were 
zero, the initial condition would be imposed upon another yf k \ not all of which 
can vanish at t = 0. The arbitrary can be restored if desired in the final solution 
by multiplying by any factor not zero. We shall consider only enough steps 
of the process to enable us to exhibit its peculiarities and to establish its general 
character. 

Terms Independent of (i. The terms of the solution independent of (i are 
defined by the differential equations 

n 

(?£>)' + «f yS - 2 a i} y$ = 0, i=l,....,n, 

the general solution of which is 

^» = i^ Cy 6<f-^ t=l,....,n, 

j=i 

where the vjjf are the constants of integration. Since the y$ are to be made 
periodic with the period 2n, and since aj 0) — af^Oraod V — 1 by hypothesis 
except when j = k, it follows that >j$? = when j df: h. Since y lk (0) = c lk for 

all | (i | sufficiently small, it follows that yf£ (0) = c lk , y$(Q) = 0,j—l, , oo . 

Therefore we must make Y\ kk = 1. The solution satisfying the conditions laid 
down is therefore found to be 

y?k = c ih , » = i, ...,«. (4i) 

Coefficients of ft. The differential equations for the coefficients of the first 
power of ft are 

n n 

W + <« - 2 a iiy $ = - afy® + 2 «, t = 1, . . . . , n. (42) 

3=1 J=l 

When the left members of these equations are set equal to zero their general 
solution is 

^ = 2»$V <a?,, - of,)t , »=1, ...., n, (43) 

.7=1 

where the q$ are the constants of integration, and the c tj are the same as before. 
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We shall find the complete solutions of (42) by the method of the variation 
of parameters. Regarding the >7J" as variables and imposing the conditions that 
(43) shall satisfy (42), we get 

2 (>$) V (a > 0> - a?,)t = - affj/S + 2 « = fflt), i = l,....,n, (44) 

i=l 3=1 

where the g$(t) are periodic in t with the period In. The determinant of the 
coefficients of the (»$)' is 

. ,,2 (a<»)-a(0))t 2 (a(0)-a(0))t 

A = \c tJ \e*= 1 ' * =e'" 1 ' * » 

which can not vanish for any finite t. Therefore the solutions of (44) for the 

(*$)' are 

(^))' = e -(«f-«f»*Aji>, (45) 

where the AjjS? are known periodic functions of t with the period 2n. 
For j df. h the solutions of (45) have the form 

nJS = c-ci-MP)* P jk (t) + B$, (46) 

where the Pj]P (t) are periodic with the period 2 7t and the B$ are arbitrary 
constants. For j = Jc we have, from (44), 

MK)' = 4B = -aJP + «a, (47) 

where 5$ is A^ after the terms — a^yfk have been omitted from the &th 
column. It is a periodic function of tf with the period 2n and has in general 
a term independent of t. Hence we may write 

where d^ is constant and the mean value of Q$(t) is zero. Then (47) becomes 

W = (4 1) -4 1, )+ QW). (48) 

In order that qftl shall be periodic we must impose the condition 

oj? = 4", (49) 

after which we get 

*% = p& + m, (so) 

where P$. is periodic with the period 2 n and J?j^ is the constant of integration. 
Substituting (46) and (50) in (43), we have as the general solution of (42) 

n n 

2$ = 2 B$c^f>-^+ Zc if P%{t), i = 1, . . . ., n. (51) 

11 
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In order that the y$ shall be periodic with the period 2 n, all the B$ must be 
zero except B$.. From y$(0) = 0, we get 

i?o> = - J- 2 Cl ,pw (0). (52) 

Therefore the solution satisfying all the conditions is 

2$ = I \c v P$ (0 - C f- qy P$> (0)1 . (53) 

It remains to be shown that the integration of the coefficients of the higher 
powers of [i can be effected in a similar manner. Suppose aj?, . . . . , a^ _1) and 
the y$, y$, . . . . , «4™ _1) satisfying the differential equations have been uniquely 
determined so that the y { $. are periodic with the period In and that y$(0) = 0, 

1 = 1, , m — 1. We shall show that the yffi can be determined so as to 

satisfy the same conditions. 

From equations (39) and (40) we find 



n 

„(0) 



WY + <4 0) y\t - 2 a* 3#° = - < yt + 2 6%» y$ 

to— 1 n 

+ X [- a?> y%~ p) + 2 Off 2/jr- p) ] . (54) 

p=l j=l 

Omitting the terms included under the sign of summation with respect to p, 
these equations are identical in form with (42) except that we now have the 
superscript (m) instead of (1). The integrations proceed as in the case treated, 
for the terms included under the summation with respect to p are all periodic 
with the period 2n and do not change the character of the gffi (t). Therefore 
ai m) and the yffi can be uniquely determined so as to satisfy the differential 
equations and be periodic in t with the period 2n, and so that at the same time 
y$°(0) =0. Hence the induction is complete and the process can be indefinitely 
continued. 



§11. Construction of the Solutions when af } — o4 0) = mod »/ — 1. 

Suppose the aj 0) are all distinct, that <4 0) and af> differ by an imaginary 
integer, and that there are no other such congruences among the aj 0) . The 
solutions associated with <4 0) , .. .., a^ are computed by the methods of §10. 
It was shown in § 6 that in this case, in general, a lf a z and the y n , y i% can be 
developed as converging series in integral powers of p. It will be assumed 
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that we are not treating one of the exceptional cases where the series proceed 
in fractional powers of jx. 

The general solution of (39) for the terms independent of (i is in this case 

y^^infc^-^, i = i,....,n, 

where the yjffl are the constants of integration. 

Imposing the conditions that the y$ shall be periodic with the period In 
and that yii ) (0) = c u , these equations become, since af — af* is an imaginary 
integer, 

¥ff = (l - iff ^) «a + irf? c iZ eH*-«?n, » = !,....,», (55) 

where 5?ff is so far arbitrary. 

Coefficients of p. It follows from (39) and (40) that the coefficients of (i 
must satisfy the differential equations 

n n 

<lfty + aPlft-2a il ift = - a Pift + 2Wtft, i=l,----,n. (56) 
The general solution of these equations when their right members are zero is 

^ = lyi$>c v e^-^\ •=!, ....,n. (57) 

Considering the coefficients ^ as functions of t and imposing the conditions 
that (56) shall be satisfied, we get 

X «)' '*,**•-«"* = - afVff + 2 flfljfl, • = 1, . . . ., n. 

Substituting the values of y$ from (55) and solving, it is found that 

«)' = - < (i - nff J) + iff A£>(0 + b® (t), 



frff)' = - «?' iff + ^ Aff (0 + A°l (0. 



(58) 



where the AjJ 1 and Z)^' are periodic functions of t with the period 2it depending 
upon the 6$ and ^f-f)', In the first two equations the undetermined constants 
a§° and Vff enter only as they are exhibited explicitly. 

Equations (58) are to be integrated and the results substituted in (57). In 
order that the y$ shall be periodic we must impose the conditions 
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-af> (i-^) + *ff*ff + « = <>, 

-^ntl + n$m + d$ = o, I" (59) 

where 5^, J^, eZff, d^ are the constant terms of Ag } , Aff, D$ and J9ff respect- 
ively, and where the B$ are the constants of integration obtained with the last 
n — 2 equations. Eliminating yj$ from the first two equations, we get 

af - (iff + d® + J dff) a?> + (6g> dg> - Off ^') = 0. (60) 

There are two cases, according as the discriminant of this quadratic is not zero 
or is zero. In the first case, which may be regarded as the general case, the two 
roots for a^ are distinct, corresponding to distinct values of \ and b % given in 
(30) in the existence proof. It was shown there that in this case the solutions 
proceed according to integral powers of p. In the second case, corresponding 
to 5j = h z , the character of the solutions depends upon the coefficients of terms 
of higher degree, and they may proceed according to powers of (i or ± V (i. 
We shall assume that the discriminant is distinct from zero and proceed to the 
construction of the solutions. 

It will be shown that after one of the two pairs of values of a| u and riP 
satisfying (59) is chosen, the solution is uniquely determined except for the 
arbitrary constant factor which may be introduced at the end. Integrating (58), 
substituting the results in (57), and determining the arbitrary constants so that 
the solution shall be periodic, and imposing the condition that ^u(0) = 0, we find 

»ff = B$[- *** + % *M>*] + £ [c w P$ (0 - 2U c a Pj« (0)] , (61) 

» = 1, ...., n, 

where B$ is an undetermined constant and the P$ are entirely known periodic 
functions of t having the period 2?i. 

Coefficients of fi z . The coefficients of [j? are defined by 

(0ffy + «M - 2 % «ff = - <M - <M + 2 [« + «]. (62) 

i=l, ....,«. 

The general solutions of these equations when the right members are put equal 
to zero are the same as (57) except that the superscripts are (2) instead of (1). 



(63) 



(64) 
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Varying the ^jf, we find for the equations corresponding to (58) 

«y = _ of > (l - ni J) + of> m ~ + B& 4B> (0 +. Z)j? (0, 

( n $y = - of> «ff + < B$ + B® Ag> (") + Z><f (0, 
[r$y = e -(-f>-i»»t Af (0, y = 3, . . . . , n. 

The undetermined constants af } and i^J' are written explicitly in the first two 
equations, and it is to be noted that A$(i) and A$(t) are precisely the same 
functions of £ as those which appear in (58). 

In order that when we integrate (63) and substitute the results in the 
equations corresponding to (57) we shall have a periodic solution, we must 
impose the conditions 

U - i® of" + (iff - af>) B$ + dff = 0, 

l?jf = 0, y = 3,....,n 

6ff, 6ffi 4ffii 4i bein g the constant terms of Aft, Ag>, £g> and Z$> respectively. 
The first two equations are linear in af and B$ and determine these quantities 
uniquely, provided their determinant is not zero. The determinant is 

= ap + i!ff(6ff + «ff£!)-aff. 

Eliminating a^ and jygf by means of (60) and (59), we get 

A = ± V~Z>, (65) 

where Z> is the discriminant of (60) and is by hypothesis distinct from zero. 
Therefore the solution of (64) for af } and B$ is unique. The sign before */ D 
depends upon which of the two roots of (60) is used. Now integrating (63), 
substituting the results in the equations corresponding to (57), and imposing 
the condition that yfl (0) = 0, it follows that 

Vf = Bg [- 1% + c^-f »<] + ,5 \c i} Pff (t) - C ^c a Pjf (0)] , (66) 

* = 1, ....,«, 

where B$ is as yet an undetermined constant. 





5(i) + a u>^L 


mi > 


7,(1) „a) 

°21 a l 
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Now consider the general step in the construction of the solution. We 
shall have 

j—i * 
where the F} v) are known periodic functions of t. When we put the right 
members equal to zero, the general solutions are the same as (57) except that 
the superscripts are (v). By varying the constants of integration, we get equa- 
tions (63) except that the superscripts are (i>) and (v — 1) instead of (2) and (1) 
respectively. The conditions for periodicity are similar to (64) and have the 
same determinant. Consequently at this step c4"> and I?£[ -1) are uniquely deter- 
mined. Hence it is evident that the process can be carried as far as is desired. 

For congruences of higher order analogous methods are applicable, and in 
the exceptional cases to this treatment the existence proof furnishes a sure guide 
for the construction of the solutions. 

§ 1 2. Construction of the Solutions when c4 0) = af>. 

For simplicity, suppose c4 0) = a£ 0) and that there are no other equalities 
among the af\ and that no two of them differ by an imaginary integer. The 
only variations from the method of § 10 are in the construction of the solutions 
associated with af. It was shown in § 7 that, except in special cases, the 
solutions are in powers of ± V p. We shall assume that we are not treating 
one of the special cases. Hence we have 

on = of + aiV + a fV + •••■, 
a 2 = of — a$V + «FV — •••-, 

s^yff + j^V + ygV + ■••-, 
jfa^ff-ygy + ygV- 

Terms Independent of p. The terms independent of (i are defined by 
(y\f)' + of yg> -£%yf = 0, i = 1, . . .,, n. 

The general solution of these equations is 

n 

t/S = *ff % + vS (% + tea) + s nJS ««•-*>*. 

.7=3 

Imposing the conditions that the yJJ* shall be periodic with the period 2n and 
that 2/g' (°) = c n > we g et 

y|S> = c u , i = 1, ....,n. (68) 



(67) 



} (71) 
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Coefficients of {&. The coefficients of [i* are defined by 

n 

(M?)' + < ) Mi ) -2«^l ) = -a?>^> = -ai 1 > Cfl , t = l, ....,«. (69) 

*=i 

The general solution of these equations when the right members are zero is 

yi\ > =^c a + ^(c i2 + tc il ) + X^c ij e^-^ } t '=i, ....,„. (70) 

By the variation of parameters we get 

n 

Solving these equations for the (rin)'> we nn( ^ 

WPy = 0, y=2, ...., n. 

Consequently 

Substituting these values of 57JP in (70), we get 

n 

y®=lB®-«? ) tlc a + B«nCi Z + tc il ) + 2B$c ij e^-^\ i=l,....,n. (72) 

J=3 

Imposing the conditions that the y$ shall be periodic with the period 2 71 and 
that y$ (0) = 0, we have 

5j}) = 0, y = 3, ...., n, L (73) 

Then equations (72) become 

Vii = (-? c « + O"^ » = 1, • • • ., n, (74) 

where a^ remains as yet undetermined. 

Coefficients of [i. The coefficients of /* satisfy 

Gfffy + ofyff-S^y^s-afyff-aPj^ + lflg)^, i = l,....,n. (75) 
.7=1 >=i 

The solution of these equations when the right members are zero is of the same 

form as (70), and we find, by varying the constants, 

M?y 0a + (nff y (<fc + * %) + 2 W <fc e<»f -P> * = -< > *,<»> - of> yg> + I flg> $>. 

j=3 j=l 
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Solving these equations, we get 

(uff y = - «f + ( J + < 2 + « A n (0 + *% w, 

w?y = u -<xi i)3 - Aino, 

(„jf )' = e -<"f -«?»< Aj?> (t), / = 3, . . . . , n, 

where Ajf (<) and Z)ff (£) are periodic functions of t. The first of these equations 
gives rise to integrals of the type 



(76) 



°*f 



. sin ., ,. a* . cos .. . (tj sin ., 

cos J j sin J p cos J 



The second equation gives rise to the corresponding integral 

When we substitute these results in the equations corresponding to (70), we get 
for these terms 

sin 



m +*«>«. =f[=FO±«£# 



cos ., 



sin ' 



! ] + * 



cos 



jt- 



Consequently the terms of the type t ■ jt vanish when the results of the inte 



sin 



gration of (76) are substituted back in the equations corresponding to (70). 
Hence, we get at this step 



»ff = 38* + 3ff (fl. + tc a ) + 2 -B^rfF-P'' 



JF=3 



+ 



[(- < + g «i« a + <2ff) * + i (af 2 + 6g>) * 2 + J 3 *? (*)] On 



+ [- (of + iff) < + /fi> (<)] c i2 + X c tj Pjf (t), 

j=3 



(77) 



where JJf and df} are the constant terms in Aff and 2){f, and where the Pjf(t) 
are periodic functions of t. In order that the y$ shall be periodic and 
yf) (o) = 0, we must impose the conditions 






■fiff = o, 



3, .. 



-Bff c n + 5g> c 12 + 2 Pjf (0) oy = 0. 
^=i 



(78) 
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The only undetermined constant remaining in (77) is af, and the y$ now have 
the form 

¥ff = (^-°fc^aT + ^(t), i = l,....,n, (79) 

\ c ii ' 

where the <I>|f (£) are known periodic functions of t. 

The constant a^ has a double determination except in the special case 
where b$ = 0. We shall suppose b$ ^ °> f° r if i* i s zero we have one of the 
special cases excluded above. It will be shown that the work becomes unique 
after one of the two possible values of aP is chosen. 

The coefficients of n§ are determined by 

(j#y + oj» s# - i % kjp = - of> »§» - «f sff - < »ff + 2 eg' «ff • 

The equations corresponding to (76) are 

w?y = - op> + 2 (g + «p «f ) + < ^ff (0 + aff w, 

(Uffy= U -2aJ»of>- Ag>(0, 

(,$> )' = e -<<f >-«f»t a^) } y = 3> . . . . ? n . 

In order that the final solutions at this step shall be periodic, we must impose 

the condition 

2a < 1 1) af + 5g ) = 0, 

which determines af uniquely since af } dp 0. The other constants are all 
uniquely determined by the periodicity condition and the initial condition except 
af } , which is fixed by the periodicity condition at the next step. Similarly, 
another solution is obtained using the other determination of aj 1 '. The solutions 
associated with c4 0) , . . . . , a^ are obtained by the method of § 10. 

The chief types of cases have been treated, and the exceptions to them are 
developed similarly, according to the forms indicated by the existence proofs. 

§ 13. Solutions when the $% do not All Reduce to Constants for (i = 0. 

Heretofore we have considered linear differential equations whose coefficients 
become constants when ^ = 0. We shall now waive this restriction, but we 
shall suppose the solutions are known for (i = 0. Let the equations under 
consideration be 

xi = I f 5 0$W~U, • = 1, . . . ., n, (80) 

12 
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where the 8$ are periodic functions of t with the period 2 n, and where not all 
the 6$ are constants. For fi = these equations reduce to 

n 

(af)' = 2 0§>af, » = 1, ....,«. (8*1) 

The solutions of these equations in general have the form 

n 

af = 2 Jf e">° )(t) g$ (t), i=l, ...., n, (8 2) 

where the y$ are periodic functions of t with the period 2th. 

Suppose the aj 0) and y<$ are fully known. We desire the solutions of (80). 
It follows from the general results of § 3 that there is at least one solution 
of the form 

Xi = e at y i , i = l, , n, (83) 

where the y i are periodic with the period 2n. Now equations (80) can be 
integrated as power series in ft, reducing to xf> for ft = 0. We form n solutions, 
x a , . . . ., x in} defined by the initial conditions x l £> = 0, x<® = 1, k = 1, . . . ., n. 
Then any solution can be expressed linearly and homogeneously in terms of 
these n solutions. Hence 

n 

x % = 2 AjXq, i = 1, . . . . , n. (84) 

3=1 

Transforming (80) by (83), we get 



yi + ay i = 2 \ 5 0$V|&> 

j=l L k=0 J 



1 = 1, 



Consequently the conditions that the y i shall be periodic with the period 2n 
are y^n) — «a(0) = 0, i= 1, . . . . , n ; or, because of (83), (84) and the initial 
conditions imposed on the x ijf 

2A j [e- 2 *"xt l (27t) — x tj (0)'] = 0, » = 1, n. (85) 

Now, since x v = 2 x[f [i k , the determinant of the coefficients of the J.,- set 

Jc=0 

equal to zero is 

A = | e-^ asg? ( 2 *) — a»g» (0) + «-*• I a$> ( 2 *) p* | = 0. (8 ft) 

fc=l 

This is an equation for the determination of a. For ^ = the solutions are 
a = <xj 0> + v V — 1, where v is any integer. Consequently, if ai 0) is a simple 
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root of (86) and if no two of the ocj 0) differ by an imaginary integer, then (86) 
can be solved for (a — a£ 0) ) as a converging power series in (i. The results sub- 
stituted in (85) give the ratios of the A } as power series in fi for k = 1, . . . ., n. 
When these results are substituted in (84), we have the solutions expanded as 
converging power series in (i, and they have the form (83), where a and the y t 
are power series in p, and where the latter are periodic with the period 2 71. 

There are other cases where for fi = the roots of (86) satisfy different con- 
ditions. The d-iscussion is parallel to that where for (i = the 6 V are constants. 
The essentials are that the differential equations shall be expansible as power 
series in a parameter (i, and that for (i = the solutions shall be known. The 
process is fundamentally one of analytic continuation of the solutions with 
respect to the parameter (i, and it can be repeated and continued from one value 
of ft to any other, provided the series do not pass through a singularity in the 
interval. 

Non- Homogeneous Linear Differential Equations.* 

§ 1 4. Case where the Bight Members Are Periodic with the Period 2 n 

and the a,j Are Distinct. 

The problems of celestial mechanics generally lead to sets of differential 
equations having right members which are functions of the independent variable 
alone. The character of the terms in the solutions introduced by these right 
members depends not only upon the properties of the latter but also upon the a 5 . 
We shall treat the most useful cases. 

Suppose 

n 

xL — 2 6 tJ x J = g t (t), * = 1, , n, (87) 

j=i 

where the d^ and g t are finite and continuous periodic functions of t having the 

period 2 71. For the left members set equal to zero the general solutions are in 

this case 

n 

Xi — 2 vuftftj, i — 1, , n, (88) 

where the y^ are periodic with the period lit and the yi } are the constants of 
integration. 

* For a different treatment of this subject see paper by W. D. MaoMillan, Transactions of American 
Mathematical Society, Vol. XI, No. 1, p. 85. 
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By the method of the variation of parameters we find 

n 

2 vje"' 1 y i} - g t (t), i = 1, . . . . , n. (89) 

3=1 

The determinant of the coefficients of the yjj is the determinant of the funda- 
mental set of solutions, and can vanish for no value of t for which the 6 tJ are 
regular (§ 2). By hypothesis the $ v are regular for all finite values of t. There- 
fore this determinant can not vanish, and is 

where A is the determinant of the y tJ . Consequently the solutions of (89) 
have the form 

>li = er*£, (90) 

where A,- and A are periodic functions of t with the period 2n, and moreover 
A vanishes for no finite value of t. Hence A,/A can be expanded as Fourier 
series of the form 

~1 = a ( j) + 5 [<> cos rat + V£ sin mi]. 

^ m=l 

Since the A, and A are power series in (i and in general A does not vanish for 
fi = 0, this result can also be arranged as a power series in (i whose coefficients 
are periodic with the period 27t, but it is more convenient here to regard it 
simply as a Fourier series. 

If <xf + w* 2 =£ 0, j = 1, . . . . , n, m = 0, , eo , we have 

e a * t ^dt = S- e i* + ei* 2 } m 9 . — ^- cos mt 

A Of m=iL af+rn? 

a] + m 2 J " 

or, 

y lj = B j + e-^P i {t) } (91) 

where the Pj(t) are periodic with the period 2n and the JS, are constants of 
integration. Substituting in (88), we get 

n n 

x t = 2 ^«**y„ + 2 y v P,(«), » = 1, . . . . , n, (92) 

^■=i 3=1 

as the general solutions of (87), the Pj(t) being periodic with the period 2 n. 
The terms y^Pj(t) are the ones due to the presence of the g^t). 
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Now let us suppose that a k = v y/ — 1, where v is an integer. Then 

the term 

f e -vV=it [ av cog v f _^_ i^ s i n v {\ dt 

becomes 

1 (a, — b v V — 1) £ + — (a„ V — 1 — 6„) (cos 2 vt — V — 1 sin 2 vt). 

Hence in this case we get 



x i = ^B } ^y id + l l y v P j {t) + \{a v —hW--i-)te^ l y u , i = l,....,n. (93) 



Hence, »/ £Ae a$ are distinct and none of them congruent to zero mod V — 1, and if 
the g i (t) are periodic with the period 2n, then the particular integrals are also 
periodic with the period 2n; but if one of the characteristic exponents is congruent 
to zero mod */ — 1, then the particular integrals in general will contain, in addition 
to periodic terms, the corresponding parts of the complementary function multiplied 
by a constant times t. 

§15. Case where the Right Members Are Periodic Terms Multiplied by an 
Exponential, and the a,j Are Distinct. 

We now suppose the g t (t) have the form 

ft (*) = «*/«(*), (94) 

where the f{t) are periodic with the period 2tc. When A, is a pure imaginary, 
•/I = I V — l, the gi(i) have the form 

CO 

g t (t) = 2 [a k cos (k + I) t + b k sin (Jc + I) t], 

which appears often in applications to celestial mechanics. 
If we transform the differential equations 

n 

x t — 2$ v x J = J t f t (t) (95) 

by 

x t = e H z f , 
we obtain 

4 + Jlz, — 2 d v Zi =f t (t), i — \, , n, (96) 

3=1 

which have the same character as the equations treated in § 14. If the a 3 are 
distinct, then the characteristic exponents ctj — X belonging to (96) are also 
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distinct. If no a 5 — X is congruent to zero mod V — 1, then the general 
solutions of (96) are 

where the Qi(t) are periodic with the period 2n. 

Therefore in this case the general solutions of (95) are 

xt = 2 Bjtfityy + e* Q t (t), i = l, ...., n. (97) 

But if X is congruent to one of the a jf say <x„, mod V — 1, then the z { 
have the form 

^=1 .7=1 

and therefore 

x, = 2V t y„ + « M 2y»P J (0 + i («,- 5, V=l) <«**?„. (98) 

Therefore, i/" $e # { (£) are e At ta'mes periodic functions and if no a,- — X is congruent 
to zero mod s/ — 1, then the particular solution is e Kt times a periodic function ; but 
if a„ — X is congruent to zero mod */ — 1, then the particular solution is e Kt times 
a periodic function plus a constant times t e** l y iv . 

§16. Case where the Right Members Are Periodic and a 2 = a 1 . 
In case there are no equalities among the c^ except a 2 = a x , the solutions of 

n 

sci — 2 6 v Xj = 0, i = 1, , n, 

in general have the form 

as, = y !l e^ t y il + %«"«* (y i2 + ty a ) + £ »&<*** jfo, i = 1, . . . . , n. (99) 

j=s 

For the non-homogeneous equations 

n 

xi — X6 iJ x j = g i (t) 
j-\ 

we find, by the variation of parameters, 

e^ynn'i + e ttl * (y« + ty a )n 2 + I «**y ^ = g t (t), » = 1, . . . . , n. 

J=3 
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n. 



Solving these equations for the v\\ , we get 

A >7i = | 9i (t), (y i2 + t y a ), Via, ••■■, Vin\ e~ axt , 
A*7a = | Viu & (t), Via, •••, Vm\ e - " 1 ' , 

A>7; = | ttu , {Vi* + 'y*i)i 3te, — , y in | e~ a >' , y = 3, . . 

where A is the determinant | y tj \ . The expansions of these determinants have 

the form 

n[ = e-^P x (t) — e-^H P z (t), 

V ' % -e-^P,{t), I- (100) 

n} = e->*P J {t), j = Z,....,n, 

where the Pi(i), . . . ., P n (t) are periodic with the period 2n. 



If no <x,j is congruent to zero mod *</ — 1, then we find 

m = e~^B 1 {t) — e-^HBt (t) + B, , 
y, =e~^B z {t)+B 2 , 
1,=e-**Rj(t)+Bj, j = 3, ....,«, 



(101) 



where B 1} . . . . , B n are periodic with the period 2n. Substituting (101) in (99), 
we get 

n n 

x, = B x e^y ix + B z e^(y iz + ty a ) + 2 B^y v + 2 22, y„. (102) 

.7=3 ?=1 

The terms introduced into the solutions by the right members of the differential 

n 

equations are 2 Bji/y, which are periodic. Therefore it follows that equalities 

among the characteristic exponents, without congruences to zero mod \f — 1, do not 
introduce non-periodic terms into this part of the solutions. 



The case where one a,, j = 3, . . . ., n, is congruent to zero mod V — 1, 
is a combination of the second part of §14 and this case; and that where a 2 = a 
is congruent to zero mod V — 1 does not differ from that where a z = a t = 0. 

Consequently we consider the case a z = a x = 0, which frequently arises in 
celestial mechanics. Then the equations corresponding to (100) become 

y I >=P 1 {t)-tP z (t), 

Vz=P Z (t)> [(103) 

*6 = e-* Pj(t), / = 3, ...., n, 
where 

P } = 2 |"fl#> cos Jet + btf sin kq, / = 1, ...., n. 

k=0~ 
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Hence we find 

m = a$H — \a$>t* - tB z (t) + R^t) + B u 

yi z = a<*H + R z (t) + B s , 1(104) 

n^e-^R^ + B,, y = 3, ....,«, 

where the Rj{t) are periodic with the period In. Substituting these values 
in (99), we get 

«!i = B lVa + B,(y iz + ty a ) + 1 S^'y* + (a$H + i«M y a 

+ a$Hy is + ZR j (t)y ij . (105) 

Hence, M?#era the g t (t) are periodic with the period 2 n, and two of the a,j are not 
only equal but also equal to zero, then the particular integral involves not only t but, 
in general, also t z outside of the trigonometric symbols. It can be shown similarly 
that when k of the a,- are equal to zero, then in general the solutions are poly- 
nomials in t of degree h whose coefficients are periodic functions of t. 

The University of Chicago, May 25, 1910. 



